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ABSTRACT 

A new  borate  material,  BiB306  (BiBO),  has  been  reported  to  have  extremely  high  nonlinear  optical  (NLO)  coefficient.  In 
our  investigations  we  obtained  both  glass  and  single  crystals  of  BiBO.  Owing  to  its  composition,  the  material  is  very 
viscous  when  molten.  This  factor  causes  serious  problems  both  during  synthesis  and  crystallization  of  BiBO.  In  this  paper 
only  properties  of  BiBO  glass  are  reported  We  found  that  for  the  increase  of  the  long-range  ordering  the  photoinduced 
optical  second  harmonic  generation  (SHG)  shows  a maximum,  which  is  shifted  towards  the  lower  photoinducing  power 
with  the  increase  of  the  temperature.  Maximal  value  of  SHG  was  about  0.081  pm/V  at  1.06pm  radiation  from  YAG:Nd 
laser  (W=30MW,  r=35ps). 

keywords:  borate  glass,  bismuth  borate,  nonlinear  optical  properties. 

1.  INTRODUCTION 

In  recent  years  borates  have  been  intensively  investigated  due  to  their  excellent  nonlinear  optical  and  laser  properties. 
Among  the  newly  synthesized  borate  single  crystals,  it  is  necessary  to  point  out  the  following:  BiBO,1  YCOB,2  KAB,3  and 
CLBO.4  Crystalline  BiBO  was  reported  by  Liebertz  et  al.1  to  have  extremely  high  NLO  coefficient  for  SHG  of  3.2  pm/V. 
This  value  is  larger  than  that  of  widely  used  NLO  crystals  such  as  p-BaB204  (BBO)5  or  KTiOP04  (KTP).6 

In  our  investigations  both  BiBO  crystals  and  glass  were  produced.  In  this  paper  we  report  only  on  photoinduced 
optical  effects  in  BiBO  glass,  the  NLO  properties  of  BiBO  single  crystals  will  be  published  elsewhere.  Molten  borates  are 
characterized  by  very  large  viscosity  and  its  steep  temperature  dependence.7  Some  molten  borates  have  viscosity 
approaching  103  poise  near  their  temperature  of  crystallization.  The  segregation  of  melt  constituents  on  the  melt-crystal 
interface  is  very  slow  and  even  the  slightest  changes  of  the  melt  temperature  can  deteriorate  the  growth  conditions,  what  in 
consequence  gives  crystals  with  poor  quality.  Additionally  such  melts  show  strong  tendency  to  glass  formation,  so  it  is  very 
difficult  to  initiate  their  crystallization.  On  the  other  hand  this  property  of  molten  borates  makes  obtaining  them  in  form  of 
glass  relatively  easy.  Decreasing  the  temperature  of  molten  borates  below  the  melting  point  causes  rapid  increase  of  melt 
viscosity,  what  reduces  the  atomic  mobility  in  the  melt.  The  nucleation  of  crystallization  centres  in  such  melts  is  very  slow 
and  as  a consequence  glass  formation  is  observed. 

The  starting  material  was  synthesized  from  stoichiometric  amounts  of  B203  and  Bi203  in  a platinum  crucible.  To 
obtain  material  with  known  composition  boron  oxide  was  molten  first  to  remove  water  which  is  easily  absorbed  by  this 
compound.  In  next  step,  after  finding  the  mass  of  molten  B203  appropriate  amount  of  Bi203  was  added  to  the  crucible. 
Owing  to  high  viscosity  of  molten  BiBO  it  was  very  difficult  to  obtain  homogeneous  melt  which  had  to  be  stirred  with  the 
use  of  a platinum  mixer  at  850°C.  The  process  lasted  several  days  until  the  melt  became  totally  transparent.  BiBO  glass  was 
produced  by  rapid  cooling  of  the  melt  to  temperatures  below  600°C,  then  rods  of  different  diameters  were  formed  from  the 
viscous  melt. 

For  investigations  we  have  chosen  the  specimens  with  different  degree  of  crystallinity.  The  degree  of  crystallinity 
was  evaluated  as  a ratio  of  intensity  of  the  first  crystalline  maximum  to  the  total  amorphous-like  background.  The 
corresponding  parameter  was  introduced  by  direct  ratio  of  the  first  X-ray  maximum  to  the  total  amorphous-like  background. 

2.  MEASUREMENT  PROCEDURE 

The  output  optical  SHG  light  intensities  were  detected  using  FEU-79  (FEU-39)  photomultipliers.  The  measurements  were 
carried  out  in  the  single-pulse  regime,  with  a pulse  frequency  repetition  of  12  Hz  with  the  tunable  pulse  duration  within 
1.0  ps.  From  expression  T«l-pilp,  where  p-is  the  two-photon  absorption  (TPA)  coefficient,  1-  is  the  thickness  of  specimen 
(about  0.5  mm),  T-intensity  dependent  transparency,  we  have  determined  the  TPA  coefficient.  The  nonlinear  TPA  was 
extracted  from  the  intensity-dependent  transparency  measurements.  Simultaneously  we  have  subtracted  an  influence  of 
linear  absorption  (220  cm"1,  k=530  nm)  and  of  Fresnel  reflection  0.04.  The  averaging  statistics  over  the  sample  surface  was 
performed  in  order  to  average  a space  nonuniformity  in  the  optical  constant  distribution. 
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During  the  evaluations  of  time  delayed  nonlinear  optical  response  we  performed  measurements  of  the  light 
intensities  at  © and  2©  laser  frequencies  with  the  time  step  of  about  0.06  ps.  This  was  ensured  using  the  electronic  boxcar 
(GENETIC  12  M)  used  in  the  time-synchronized  pump-robe  conditions.  The  calculations  of  the  nonlinear  optical 
susceptibility  were  done  from  the  expression: 
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Here  Ro  is  a radius  of  the  pumping  beam  which  possesses  gaussian-like  form;  n(2m)  and  n(<o)  - refractive  indices  for  the 
pumping  and  PISHG  doubled  frequencies  during  the  photo-pumping;  Xijk  - components  of  the  second  order  nonlinear 
optical  susceptibility  determined  for  different  angles  of  incident  light. 

To  exclude  an  influence  of  the  hyper-Raman  scattering,  we  have  performed  the  additional  investigations  of  the 
scatterred  light  up  to  3000  cm'1,  starting  from  the  wavelength  light  of  530  nm.  Three  hyper-Raman  maxima  between  1200 
cm'1  and  1800  cm'1  were  observed  with  the  intensities  at  least  13  times  weaker  than  the  probing  phototransparent  signal.  In 
the  case  of  the  Raman  and  hyper-Raman  measurements  we  have  used  argon  laser  at  522  nm  Ar+  laser  line  and  Spex 
Triplemate  spectrometer.  A position-sensitive  photomultiplier  calibrated  using  helium-neon  gaseous  mixture  discharge  was 
used  to  detect  the  scattered  light  quanta. 

Photoinduced  changes  of  the  investigated  specimens  were  done  by  light  of  the  Q-switched  Ti-saphire  laser.  The 
pulses  of  this  laser  were  multiplied  by  frequency  and  their  power  was  varied  within  the  0.1  - 12  MW,  the  corresponding 
pulse  duration  was  changed  within  the  2-  65  ps.  As  a probing  second  harmonic  generation  source  the  pulses  of  the  YAG-Nd 
laser  with  wavelength  1.06  pm  and  pulse  duration  2 ps  were  used.  The  output  second  harmonic  generation  (4=0.53  pm)  and 
pumping  (4=1.06  pm)  signals  were  spectrally  separated  by  grating  monochromator  SPM-3.  The  detection  of  the 
corresponding  signals  was  performed  by  450  ps  boxcar  connected  with  the  PDP  computer. 

From  the  Fig.  1 one  can  clearly  see  that  increasing  photoinducing  power  E and  disordering  (otherwise  connected 
with  degree  of  crystallinity  dc)  stimulate  the  enhanced  value  of  the  SHG.  This  indicates  an  additional  enhancement  of  the 
output  SHG  signal  due  to  the  disordered  non-centrosymmetric  voids. 
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Fig.  1.  Dependence  of  the  SHG  versus  the  photoinducing  power  E and  degree  of  crystallinity  dc. 
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Special  interest  causes  influence  of  scalar  hydrostatic  pressure  effectively  varying  the  inter-atomic  bonds.  From 
Fig.  2 one  can  see  that  increasing  hydrostatic  pressure  suppresses  the  output  TPA  signal.  At  the  same  time  the  decreasing 
temperature  favours  increasing  TPA.  These  values  are  described  by  the  fourth  rank  tensors  and  in  this  case  dominant 
contribution  gives  randomly  oriented  vacancies.  This  feet  is  very  important  for  possibility  of  operation  by  the  output  TPA  in 
the  partially  disordered  materials  due  to  varying  dc. 


3Q0 


Fig.  2.  Temperature  dependence  of  the  TPA  versus  the  hydrostatic  pressure  and  temperature.  All  the  measurements  were 
done  for  the  dc  = 12  %. 


Fig.  3.  Dependence  of  the  density  of  states  of  BiBO  glass  on  energy  for  the  different  degree  of  crystallinity:  • -5%;  □ - 12 
%;  0 -24%;  + - 38%. 

From  Fig.  3 one  can  see  additional  redistribution  of  the  density  of  states  for  the  amorphous-like  glasses  of  different 
ordering  obtained  from  X-ray  photoelectron  spectroscopy.  One  can  clearly  see  that  enhanced  disorder  substantially 
decreases  total  density  of  states.  In  order  to  explain  the  observed  phenomena  we  have  done  theoretical  simulations  using  the 
molecular  dynamics  simulation  and  solid  state  calculations.  Generally  the  method  is  similar  to  the  described  in  the  Ref.  8. 
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Calculations  of  the  band  electronic  structure  have  been  performed  taking  into  account  electron-phonon  interactions. 
The  calculations  have  been  carried  out  using  ab  initio  norm-conserving  pseudopotential  method  within  a Local  Density 
Approximation  (LDA).  The  degree  of  crystallinity  has  been  taken  into  account  using  an  extended  quasi-Brillouin  zone.  We 
have  separated  contribution  of  the  electron  and  phonon  subsystems,  effectively  contributing  to  the  nonlinear  optical 
susceptibilities.  The  phonon  modes  (within  the  BiBO  glasses)  have  been  calculated  using  a harmonic  approximation,  i.e.: 

d2T'k/dQk2+[47r2  ^h'2  Qk2]¥k  = 0 (2) 

where'f'k  is  a phonon-like  wave  function  corresponding  to  k-th  normal  coordinate  Qk  witihin  the  BiBO  crystallites, 
renormalized  by  surrounding  long-range  amorphous-like  background;  pk  denotes  a reduced  mass  of  ions  in  the  k-lh  phonon 
mode.  Second  derivatives  of  the  effective  crystalline  potential  (renormalized  by  surrounding  amorphous-like  background) 
play  a dominant  role  in  the  observed  phonon  modes. 

An  eigen-energy  of  the  k-th  phonon  mode  can  be  written  in  the  following  standard  form: 

Qk(Vk)=2Qko(Vk+l/2)  (3) 

where  Qko  = h/(fk/pk)1/2/2  is  an  eigen-energy  for  the  zero-point  quasi-phonons  and  vk  = 0,  1,2,  ...  are  phonon  quantum 
numbers  associated  with  the  following  wave  functions: 

M/k(Qk)=(2flk<y/  7r)  1/4(2Vk  /vk!)-1/2exp(-£2ko  Qk2)Hv((2Dk0)I/2Qk)  (4) 

where  Hv(x)  is  a Hermite  polynomial. 

In  order  to  model  the  vacancy-induced  changes,  it  should  be  noticed  that  the  electron-phonon  interactions  should  be 
considered,  particularly  for  the  low-frequency  phonon  modes  that  could  be  in  a frequency  resonance  with  vacancy 
resonances.  Calculations  of  the  electron-phonon  anharmonic  potential  have  been  carried  out  (this  can  be  seen  in  absorption 
spectra ) in  a nonlinear  approximation,  i.e.: 

V^rO  = e2!^  Miik"1/2  [Zms(rs-ums)|rs-unBr3  - Zm-s'(rs-um.s  )|rs-um.s.r3]  (5) 

where  and  eZ^  is  effective  ionic  mass  and  charge  for  corresponding  ions  numbered  by  m and  s,  respectively.  The 
intra-  as  well  as  inter-crystalline  interactions  have  been  taken  into  account.  Moreover,  two  effective  wavevectors  have  been 
introduced.  The  first  one,  traditionally  corresponding  to  the  ordered  crystallite  units  and  the  second  one,  quasi-wavevector 
that  is  built  on  the  effective  molecular  sphere  originating  from  the  amorphous-like  disorder.  Influence  of  the  interface  region 
(between  the  effective  BiBO  nanocrystallites  and  the  glass-like  background)  has  been  extracted  by  performing  a derivative 
procedure  that  eliminates  possible  potential  jumps  on  the  crystallite  borders.  The  u vector  is  a relative  displacement  of 
two  ions  from  their  equilibrium  positions  rs  and  rS'  . A probability  of  a one-phonon  transition  induced  by  the  phonon  of  a 
frequency  Qk  is  equal  to: 


H(Qk)=4(h/2u)'2c3H-1g-,(r1)(Eel-Qk  f ©(Qk)  (6) 

where  H is  a sum  of  the  r|  and  £,  level  widths,  Ed  is  an  electron  transition  energy,  Qk  denotes  a phonon  energy  and  g(rj)  is  a 
degeneration  degree  of  the  corresponding  electron  energy  levels.  The  parameter  ©(Qk)  is  equal  to: 

©(Qk)  = |{Lp<hdldVe.Ph(ri)|cp,'nJM>  <tp|d|£>  (ErEn+Qk)''  + 

+ Z^<Tiid|q»  «P,Tln!Ve.ph(r,)|^'nn.1>  (E^-Ep-Qk)'1^  (7) 

here  t]  and  £are  lower  and  upper  electron  energy  levels,  respectively;  cp  denotes  a virtual  electron  state,  d is  an  electric  dipole 
moment  for  a given  spectral  transition.  The  summation  is  performed  over  all  degenerated  initial  and  final  states. 
Orthogonalization  to  the  highly  delocalized  states  has  been  carried  out  by  the  orthogonalisation  procedure  of  the  Legandre 
polynoms.  The  symbol  0 denotes  averaging  over  all  occupied  phonon  states  for  phonons  with  frequency  Qk.  The  symmetric 
phonon  modes  included  in  the  electron-phonon  interaction  ( see  Eq.  5 ) are  renormalized  by  the  normal  coordinates  that 
leads  to: 


©(Qk)  =C^(rkA)C/(rvA’)  hn  GM."’(  rxA,  Qk2) 


(8) 
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where  GM’n  OA)  is  a Green  function  (y  and  y’  are  numbers  of  interacting  coordination  spheres)  defined  as: 

Gm'V)  = 2Z<p { <rj I V e-Ph(r 01  q>>  «p|d|£>  + <r,|d|(p>  <q>|V61*(ri)||>  J^-E,)"1  (9) 

The  Green  function  calculations  were  carried  out  for  a perfect  BBO  lattice  renormalized  by  glass-like  environment 
with  taking  into  account  effective  nanocrystallite  sizes  and  by  performing  the  summations  over  125  k-points  in  the 
irreducible  part  of  the  quasi-Brillouin  zone.  The  resulting  expression  is  given  below: 


Gm^'C  rA  Gk2)=IoKA/(  rxA)  K/(  r,a)(Dk2-Q2  - iS)'1  (10) 

where  the  coordinates  K/(  rkA)  were  obtained  for  a given  quasi-phonon  type  by  means  of  the  electron  states  averaging.  The 
vacancy-induced  lattice  perturbation  of  the  Green  function  has  been  done  using  a deformation  localization  that  allows  us  to 
use  a Dyson  relation: 

gaa.^'(1)=  gm."’(0)  + gm.™'(0)  u G^’O)  (1 1) 

where  GM,rr  (0)  and  G^^’fl)  is  the  Green  function  for  an  ideal  and  vacancy-perturbed  system,  respectively.  Coming  out 
from  the  equation  (11),  we  have  obtained  the  effective  electron-phonon  wavefunctions  for  calculations  of  the  intra-the- 
cluster  electrostatic  potential  and  corresponding  phonon  modes  before  and  after  applying  the  acoustical  waves.  Afterwards 
the  interaction  of  the  defects  with  appropriate  phonon  modes  has  been  taken  into  account.  The  changes  of  effective 
electrostatic  potential  contours  and  of  effective  low-frequency  vacancy  modes  under  influence  of  external  photoinducing 
field  show  that  the  applied  photoinducing  field  essentially  changes  the  electrostatic  potential  distribution.  The  appearance  of 
additional  non-centrosymmetry  in  the  electrostatic  potential  distribution  is  a main  requirement  to  enhance  of  the  second- 
order  nonlinear  optical  properties,  particularly  of  the  SHG,  and  such  an  effect  is  unambiguously  seen  in  this  case. 
Comparison  of  the  calculations  made  before  and  after  applying  of  photoinducing  field  indicates  on  essential  charge  density 
redistribution  before  illumination  and  after  the  illumination.  The  branches  corresponding  to  the  acentric  symmetry  are 
responsible  for  the  SHG  and  the  centrosymmetric  ones  for  the  TPA. 

3.  CONCLUSIONS 
In  conclusion  it  is  necessary  to  underline  that: 

Experimentally  and  theoretically  was  revealed  an  occurrence  of  the  non-centrosymmetry  due  to  the  disappearance  of  long- 
range  ordering  in  the  BiBO  glass. 

Increased  non-crystallinity  stimulates  an  enhanced  SHG. 

Two-photon  absorption  is  critically  dependent  on  the  hydrostatic  pressure. 

Band  energy  structure  calculations  confirm  essential  redistribution  of  the  valence  electronic  states  during  the  photoinducing 
illumination. 
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